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Strings, matrix models, and meanders*^ 

Y. Makeenko''* 

Institute of Theoretical and Experimental Physics, Russian Federation 
and The Niels Bohr Institute, 2100 Copenhagen, Denmark 

I briefly review the present status of bosonic strings and discretized random surfaces in D > 1 which seem to be 
in a polymer rather than stringy phase. As an explicit example of what happens, I consider the Kazakov-Migdal 
model with a logarithmic potential which is exactly solvable for any D (at large D for an arbitrary potential). 
I discuss also the meander problem and report some new results on its representation via matrix models and the 
relation to the Kazakov-Migdal model. A supersymmetric matrix model is especially useful for describing the 
principal meanders. 



1. Introduction 

I begin this talk with a brief review of the 
present status of the bosonic Polyakov string and 
discretized random surfaces in I? > 1 dimen- 
sional embedding space. As an explicit example 
of what happens, I consider then the Kazakov- 
Migdal model with a logarithmic potential which 
is exactly solvable for any D (at large D for an 
arbitrary potential). I discuss at the end the chal- 
lenging meander problem which is more compli- 
cated than the ones solved before by means of 
the matrix-model technique but maybe is simpler 
than the large- iV QCD. 

2. Bosonic string in 13 > 1 

2.1. The D = 1 barrier 

The D — 1 barrier is associated with the KPZ 
(Knizhnik-Polyakov-Zamolodchikov) formula [Q 
(which was in fact known ||,| before KPZ) 



7str 



D-1- - -D)(25- D) 
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(2.1) 



for the critical index of string susceptibility of the 
bosonic Polyakov string in a Z3-dimensional em- 
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bedding space. Alternatively, it describes two- 
dimensional quantum gravity interacting with 
conformal matter of the central charge c = D. 

The right-hand-side (RHS) of Eq. (|]|) is well- 
defined for D < 1, where it is associated with 
topological theories of gravity which can be de- 
scribed also by (multi-) matrix models. The RHS 
becomes complex for 1 < D < 25 which is physi- 
cally unacceptable. There are two alternatives of 
how to interpret this fact: 

• KPZ- approach is not applicable for D > 1 
(say other structures of conformal theory 
may become relevant). 

• KPZ-approach is correct but the interpreta- 
tion of the result is not straightforward (say 
the theory no longer describes a string). 

I discuss below that the second alternative re- 
alizes: the theory is in a branched polymer rather 
than stringy phase. 

2.2. The DDSW-mechanism 

The mechanism which governs the Polyakov 
string {— discretized random surfaces) in D > 1 
was discovered by Das, Dhar, Sengupta and Wa- 
dia [Q. They considered the curvature matrix 
models, where the propagator is modified as 



('J'^'J'fe;) Gauss 



(2.2) 



with the Hermitean N x N matrix A describing 
an external field. 



The modification (2.2) of the propagator 
results in the following modification of the 
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quadratic part of the potential 



tr$^ 



trA^A^. 



(2.3) 



These models are solved as —> cxd for a logarith- 
mic (Penner) potential in and for an arbitrary 
potential in 

The partition function in the external field A 
acquires the extra multiplier 



n 



(2.4) 



where is the coordination number of a given 
vertex a (which describes intrinsic curvature). 
Representing the matrix <I> in the form $ = ri^Ai7 
with diagonal A and integrating over the uni- 
tary matrix f2, one gets the action which involves 
interaction terms in the form of the products 
tr (J)*^ • • • tr with A-dependent couplings. 

The simplest such a modification of the action 
reads H 



S = 



ftr<i>2 



^tr$4 



5'tr$^tr$^ 



(2.5) 



The coupling constant g' describes a new kind 
of interaction due to touching of surfaces. These 
touching diagrams do not vanish as —> oo be- 
cause of the Weingarten arguments : the small- 
ness {N~^) of connected correlators is compen- 
sated by the fact that the action ^ iV^. More 
than one touching of the same surfaces is sup- 
pressed as 7V~^. 

As was shown by DDSW, the critical index 7str 
can be expressed at some value of g' via that 
(7c<i) without the touching interaction by the 
formula 



7str 



7c<i 



7c<i 



(2.6) 



While only 7c<i — — 1/2 is possible for the quartic 
self- interaction (2.5) with g' — 0, this formula can 
be extended to arbitrary 7c<i which is associated 



with the c < 1 KPZ-formula (2.1) 



2.3. Trees of baby universes 

Equation ( ^.6[ ) can alternatively be obtained 
by making a resummation in the sum over sur- 
faces with touching included. Typical surfaces 
which lead to the critical behavior (2.6) in > 1 
are trees constructed from closed two-dimensional 




Figure 1. Trees of 2d baby universes in the D- 
dimensional embedding space. 2d theory at eac h sur- 
face is critical with 7c<i < while Eq. (2^) with 
7str > in the D dimensions is reached by tuning the 
coupling of the touching interaction between the 2d 
surfaces. 



surfaces (baby universes) as is depicted in Fig. |l|. 
At each of 2d surfaces, the continuum limit with 
7c<i < is achieve d by tuning the cosmological 
constant while Eq. ( 2^ ) with 7str > in the D- 
dimcnsional embedding space can be reached by 
tuning the coupling of the touching interaction. 

The case of 7c<i = — 1, which is associated with 
no critical behavior at the 2d surface at all, results 
due to Eq. ( ^ ) in 7str = 1/2 — the typical mean- 
field value for pure branched polymers. 

The case of 7c<i = —1/m, which is associated 
with the standard critical behavior of 2d gravity 
(with matter), results in 7str = l/{m + 1) due 
to polymerization. It differs from the mean-field 
value 7str = 1/2 due to effects of 2d gravity. This 
reminds of how the critical index 7c<i = ~l/3 
appears for the Ising model on a random lattice. 

The formula ( |2.6| ) for the critical behavior of 
random surfaces (= Polyakov string) in the D- 
dimensional embedding space 



1) describes numerical data for c > 1 theories 
(see [H and references therein). 
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2) can be derived |9| from the Liouville grav- 
ity (by changing the gravitational dressing 
from exp (a+0) to exp (a^ (/))), 

3) is rigorously proven assuming locality, 

4) can be understood at large D (see Sub- 
sect, below). 

Thus the bosonic Polyakov string {— dis- 
cretized random surfaces) are in a branched poly- 
mer rather than stringy phase for D > 1. The 
typical surfaces are crumpled rather than smooth. 
Such a ground state is stable and has no tachionic 
excitations. This picture is expected for any ma- 
trix model describing discretized random surfaces 
in D > 1. I consider in the next Section an ex- 
plicit example of the Kazakov-Migdal model . 



3. The Kazakov-Migdal-Penner model [12] 



3.1. Three equivalent models 

A natural multi-dimensional extension of the 
matrix chain is the Kazakov-Migdal (KM) 
model [ pl| which is defined by the partition func- 
tion 

Zkm = J Y\. ^^^f n '^'^^ 



Nti 



(3.7) 



Here 



and Uxy are N x N Hermitean and uni- 



tary matrices, respectively, with x labeling lattice 
sites and {xy} labeling the link from the site x to 
a neighbor site y. 

The lattice can be one of the following: 

i) an infinite D-dimensional hypercubic lat- 
ticc 0, 

ii) a Bethe tree |l^ , 

iii) a q-simplex Q (see Fig. ||). 

These three models are equivalent as ^ oo at 
strong coupling when the coordination numbers 



A^2D 



1 



(3.8) 



coincide. In particular, the model on a triangle 
(q = 3) is equivalent to that on a one-dimensional 
chain (D = 1). The one-matrix model is associ- 
ated with q = 1 or D = 0, while the two-matrix 
model is described by g = 2 or _D = 1/2. The 




Figure 2. Lattice in the form of a 5-simplex (depicted 
for q — 5). 



gauge field Uu can be absorbed in the latter case 
by a unitary transformation of (j) so that the stan- 
dard Hermitean two-matrix model is recovered. 
The proof of equivalence is presented in ||l^ . 

The integration over the gauge field Uxy is over 
the Haar measure on SU(A^) at each link of the 
lattice. The model (3.7) obviously recovers the 
standard open matrix chain if the lattice is just a 
one-dimensional sequence of points for which the 
gauge field can be absorbed by a unitary trans- 
formation of . 



3.2. Loop equations 

A convenient way of solving the KM model is 
via the loop equations which are written for the 
one-link correlator 



xy 



(3.9) 



One gets 



=-"--^v--^^-(^(x^))<"») 

at asymptotically large v. 

The loop equation has the same form ||l^ as 
the one for the two-matrix model 



duj V{uj) 



2'Ki (i^ — Lu) 
= El, Gi,\ + c {XGuX — i?i/) 

with the potential 

V{lo) = V'{lo) - (A - l)F(w), 



(3.11) 
(3.12) 
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where F{u!) is determined by the pair correlator 
of the gauge fields 



cNtr (tiUliiiU 



(3.13) 



The contour Ci in Eq. (3.11) encircles counter- 
clockwise singularities of the function G^jx so that 
the integration over w plays the role of a projector 
picking up negative powers of v. 
The 1/A term of Eq. (|3.1lD reads 



duj V'{uj) 
27ri (A — Lo) 



E,._ 



El 



(3.14) 



which coincides with the loop equation for the 
Hermitean one-matrix model with the potential 

V\uj) = V'{lu) - F{uj) = V'{uj) - AF{lu). (3.15) 

This does not mean that the KM model reduces 
in general to the one-matrix model since V'{uj) 
may have singularities outside of the support of 
eigenvalues of the master field (t)sp whose spectral 
density p{uj) is parametrized by V{uj) as 



(3.16) 



The solution to Eq. ( [3.1l| ) consists of the fol- 
lowing steps. 

1) Given V{lo) find G^x and Ex- 

2) Calculate V"'(A) ^2YleEx at the cut. 

3) Then = AV'-(l-A)y' and F = V'-V^' 
are completely determined. 

It is evident that V — V for the one-matrix model 
(A = 0) and V = V for the two-matrix model 
(A = 1) when the solutions are known. 

3.3. Two explicit solutions 

The exact solutions of the KM model are 
known for 

• the Gaussian (quadratic) potential |l^, 

• the Penner (logarithmic) potential fl^ . 
The quadratic potential can be conveniently 

parametrized as 



^2 



{2D - l)c^ 



4>\ (3.17) 



where the parameter a/b is D-independent. The 
solution is then described by the one-matrix 
model with the potential 



b a 



which can be obtained substituting D 
Eq. (|3.17| ). The function ( |3.13| ) reads 



F(0)=c2-^ 



(3.18) 
= in 

(3.19) 



The exactly solvable logarithmic potential 
reads in the same notations 

V{4>) = -aln(6- </>) - (2D - !)(" + 1) 

x\u{a + c(t)) + [{2D-l)cb-a](j) (3.20) 

and 

a + 1 
a + c(f) 

The Gaussian solution is recovered when 

a = ab, a~6— >oo. (3.22) 

As is shown in |l^ , the quartic potential is repro- 
duced in the naive continuum limit for 13 < 4. 
Introducing new variables 



F(0) 



(3.21) 



cb- 



2c 



■I, (3 = 



a + cb 



(3.23) 



one gets 

y(0) = -aln(/3- 0) 

+ {a+l)ln{/3 + (t>)-2/3<j>. (3.24) 

The one-cut solution of the one-matrix model 
with the potential ( |3.24| ) reads ||l2|,0 



Ex = 



-V(A-a:-)(A-x+)§±^(3.25) 



2 ^'/52-A' 
where z is determined by the cubic equation 

2 



- z {0" - a - ^] + ^ 



and 



2^ ^ 2 



(3.26) 



1 ^ V(/32 - z2) (4^^ _ 1) 

x± = z ± — 3.27) 

* 2(3 2(3z ^ ^ 

are the ends of the cut. The behavior of the 
eigenvalue support and branched cuts of the log- 
arithms are depicted in Fig. ^ for various values 
of a. 
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Figure 3. Eigenvalue support of the spectral density 
(the bold line) and the branch cuts of the logarithms 
(the thin lines): a) for a > 0, b) for a —> +0, c) for 
— 1 < a < 0, d) for a — > —1 and e) for a < —1. 
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Figure 4. Phase diagra m of the KM model with 
the logarithmic potential (3.20). The bold line which 



starts at P — 0, a = — 1/2 corresponds to Eq. (3.28). 
The one-cut solution realizes for a < Uc- The critical 
lines a = ac and a = — 1 correspond to jstr = — 1/2 
and 7sir = 0, respectively, while the tricritical point 
/3 = 1/2, a — —1 is associated with a phase transition 
of the Kosterlitz-Thouless type. 



3.4. Critical behavior 

While the only continuum limit of the KM 
model with the quadratic potential is possible 
in D — 1, the logarithmic potential (3.20) re- 
veals a rich phase structure. The critical behavior 
emerges when: 

i) The spectral density ceases to be positive at 
the interval [a;_,a:+] as for the one-matrix 
models with a polynomial potential. 

ii) a;_ approaches — /3. 

iii) X- approaches x+. 

The critical behavior the type i) occurs along 
the line 



an 



(3.28) 



where a;_ = —z. At this line the discriminant of 
the cubic equation (3.26) vanishes and the one- 
cut solution is not applicable for a > ac- The 
critical behavior of the types ii) and iii) occurs 
for a = and a = — 1 respectively. These critical 
lines are depicted in Fig. |4| The critical behavior 
agrees with that for the one-matrix model with 
the cubic and Penner potentials which can be ob- 



tained from of the potential (3.24) in the limits 



^oo, {0^ - a) - P^/^ (cubic limit) (3.29) 



and 

(3 — !■ oo, a ^ \ (Penner limit), 
respectively. 



(3.30) 
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The character of the critical behavior can be 
find out by calculating the susceptibility 



f 



m Vol. da2 



InZKM 



(3.3f) 



where Vol. stands for the volume of the system = 
the number of sites of the lattice. The result in 
genus zero reads explicitly 



1 / / (/3 + x-)(/?-a:+) 
4 l^Y [fi + x+)[(5 - X-) 

2 ■ 



±1 



{l3 + x^){l3-x+) 



+Dln ■ 



(/3 + 




(/? + 


x+) 




x+) 


(/3- 





± 



•(3.32) 



where the positive sign in ± corresponds to a > 
while the minus sign should be substituted for 
a < 0. 



Having the explicit formula (3.32) for xo: we 
can find out which 7str is associated with each 



type of the critical behavior. Near the line ( 3.28| ) 
where xo is not singular and equals some value 
Xg, one gets 



Xo - Xo ^^-) for a w Uc 



(3.33) 



Since (x_ — xl) ~ {ac — aY^^ane obtains ^str = 
— 1/2 near the critical line ( ^.2g ) . 

On the contrary, xo given by Eq. ( 3.32] ) has 
logarithmic singularities for a = and a — ~1: 



Xo ^ -j^^^^iP - x+) for a w 



and 

Xo 



- log (x_ — x+f for a ■ 



-I. 



(3.34) 



(3.35) 



Near the critical lines a = and a = — 1, o ne 
gets ^str — 0. While ( 3.35| ) is positive, ( |3.34 ) is 
negative. For this reason we exclude the critical 
line a — Q from the consideration. 



3.5. Continuum limits 
3.5.1. = -1/2 

The continuum theory near a = Uc given by 
Eq. ( 3.28| ) reminds 2d gravity and can be obtained 
expanding near the edge singularity: 



4 



+ e\/A, A 



I3\ ^ 



St 



(3.36) 



where e ^ 0, A is the cosmological constant, and 
^ is the continuum momentum variable. 

The continuum spectral density is determined 
by Eq. (|3.25|) to be 



(3.37) 



and describes all continuum correlators of the 
trace of powers of the (renormalized) field 4>x at 
some point x, which is the standard set of observ- 
ables of 2d gravity. 

The critical behavior of matter is described by 
observables associated with extended objects — 
the open-loop averages 



-u. 



xy 



(3.38) 



where C'xy goes from x to y along some path an d 
the average is w.r.t. the same measure as in (3.7). 
They depend at large N only on the algebraic 
length L{Cxy) of the contour Cxy 
The double discontinuity 



C(i/,A;L) 



Bisc^I)iscxGi,x{Cxy)- (3.39) 



across the cut determines C(i^, A; 1) — the one- 
link Itzykson-Zuber correlator of the gauge fields 
which reads 



C{v,X;l) cx 



1 



(3.40) 



as e ^ 0. It is quite similar to the Gaussian one 
in the naive continuum limit at D = 1. 



The RHS of Eq. ( |3.40D describes C{v, A; L) for 
L 1/ while a nontrivial continuum limit of 
the matter correlator (3.38|) sets in for L ~ I/a/e 
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to be 

Cc{x,y;Vu) 
oc 



{x — yY + 2u(a: + y)xy + u^x^y^ 



(3.41) 



with u (X L^e. The expression ( 3.41 ) obeys the 
foUowing convolution property 



dtt^Cc{x, t; Vu)Cc{t, y; Vv) 
= Cc{x,y;Vu + Vv) 



(3.42) 



which is analogous to that ||18| for the Gaussian 
case. 

3.5.2. Jstr = 

The continuum theory near the critical line 
a « — 1 looks like 2d gravity + Id matter. The 
continuum spectral density reads 



1 



4A. 



TT 

The Itzykson-Zuber correlator 

4/32 



C(z/,A;1) 



4/32 



1 



0{e) 



(3.43) 



(3.44) 



is non-singular if f3 =/= 1/2. There is, hence, no 
unusual behavior of the matter correlators in this 
case. 

3.5.3. Tricritical point 

At the tricritical point P — 1/2 a — —I, the 
continuum system undergoes a phase transition of 
the Kosterlitz-Thouless type between the phases 
with -fstr = -1/2 and jstr = 0. 

This domain is most interesting since in the 
vicinity of the tricritical point the singular part 
of the susceptibility 



Xo 



log 



1 



1 + 3k 



2Dlog{V2S(3) (3.45) 



is Z?-dependent. Here the deviation of the tricrit- 
ical point is parametrized by 



6a = (3k + 1)(1 - K){Sf3)^ 



(3.46) 



This region is the only one where j^tr might de- 
pend on D but it was not investigated in [f2j. 



3.6. Large-D limit 

The Itzykson-Zuber integral 



I\ 



i'x,4'y] = J dUe 



cNtv<h^U''d>^U 



(3.47) 



which enters the partition functions (3.7), can 
easily be calculated as c 0: 



hil[(t>x-,(t>y\ = ciY (pxtr (fiy 



tr 

TV 



tr 

TV 



tr 



+ 0{c^). (3. 



If V{(l)x) ~ 1 as Z? ^ 00, then c - 1/L» for 
the kinetic term to be of order one and only the 
first term is left on the RHS of Eq. (|1|). The 
partition function (3.7) can be written as 



Z = 



X e 



(3.49) 



Further simplification occurs in the large-A'^ 
limit when we can replace one trace in the prod- 
uct of two traces in the exponent in ( 3.49| ) by 
the average value due to factorization. One ar- 
rives, hence, at the one-matrix model whose po- 
tential V{(j)) is determined self-consistently from 
the equation 



V{^)^V{^)-2cD{ ^ 



(3.50) 



where ( )y stands for the averaging in the one- 
matrix model with the potential V . 

For the one-cut solution of the Hermitean one- 
matrix model, one rewrites Eq. ( 3.50| ) as 



V'{X)^V'{\)-2cD 



{uj — X-){llI ~ x+) + 



(3.51) 



which is an equation for V. 

Let us identify the cosmological constant with 
gi — the coupling in front of the linear term of 
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the potential. For the susceptibihty in genus zero 
one gets 

. ^2 



Xo 



duj ■ 



Ci 



2m 



16 



(3.52) 



where / = df/dgi, while Eq. (3.51) yields 
1 



9i 



1 + cD 



(3.53) 



To obtain the critical behavior, we expand near 
X- = x'L which gives for a fc-th multicritical point 
of the one-matrix model: 



(3.54) 



Under normal circumstances when Eq. ( 3.33| ) 
holds, one gets from (3.54) 



XO - Xo - i~9i " 

so that 7str = — 1/fc since 

{9i - 9i) - {gl - 9i)- 

This is not the case, however, for 



D{x- 



x+) 



(3.55) 
(3.56) 

(3.57) 



when the denominator in Eq. (3.53) vanishes. At 
this point one has 



{x^-xl)-'^ 



(3.58) 



so that 



ist-gi) - (5j-5i)^'+'^/' - (x_-a;i)'=+i.(3.59) 

One gets, therefore, for the susceptibility 

Xo - (x_ - xir^ ^ {gl - gi)-i/('=+i) (3.60) 

which is associated with ^gt^ = l/(fc + !)■ 

The formula (3.50), which describes the reduc- 
tion of the KM model to a one-matrix model at 
large N in the larg e- £> limit, explicitly holds for 
the potential ( [3.20| ) when the exact solution is 
known at any D. As Z? — > cx) with c ^/D, 
the potential V coincides with the Penner poten- 
tial so that all the calculations can be explicitly 
done [Q. The only possible scaling behavior is 
with 7s(r = in perfect agreement with the re- 
sults of Subsect. 3.4. This seems to be a fc ^ oo 



limiting case of ^str — 1/(^ + 1) which appears 
from the critical behavior with 7str = —1/k of 
the one-matrix model with the potential V. 

Nonvanishing results for continuum correlators 
can be obtained in the large-D limit only for those 
of operators living at the same lattice site while 
the Itzykson-Zuber correlator for a contour of the 
length L is suppressed as 

C{iy,X;L) ^ ^ D-^. (3.61) 

Therefore extended correlators vanish in the 
large- Z? limit. 



4. The meander problem | 19 | 



The meander problem is known to people work- 
ing on Quantum Field Theory since the Arnold's 
question to V. Kazakov in the middle of the eight- 
ies. The problem is to calculate combinatorial 
numbers associated with the crossings of an in- 
finite river (Meander) and a closed road by 2n 
bridges .0 Neither the river nor the road inter- 
sects with itself. These principle meander num- 
bers, M„, obviously describe the number of dif- 
ferent foldings of a closed strip of 2n stamps or of 
a closed polymer chain. 

One can consider also a generalized problem 

(k) 

of the multi-component meander numbers M„ 
which are associated with k closed loops of the 
road so that M„ = Mi^\ The results of a com- 
puter enumeration of the meander numbers are 
presented in |^,|2§ up to n = 12. 

4.1. Hermitean matrix model for meanders 

Meanders can be described by the following 
Hermitean matrix model 



a=l 



Er=it'-('>^""^^»)V4.62) 



•In 



where the integration goes over the N x N Her- 
mitean matrices Wa (a = 1, . . . ,m) and 0. The 
logarithm in Eq. ( 4.62| ) leaves only one closed loop 
of the field (f). The coupling constant c is associ- 
ated with the (quartic) interaction between Wa 
and (h. 



*See 



for an introduction to the subject. 
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Expanding the generating function ( 4.62| ) in c 
and identifying the diagrams with the ones for 
the meanders, one relates the large-A^ hmit of 
^nxn{c) with the following sum over the mean- 
der numbers 



lim J- J 



NxN 



OQ 2n ^ 

11=1 fc=i 



(4.63) 



The iV — > oo limit is needed to keep only planar 
diagrams as in the meander problem. 

The RHS of Eq. ( |l6^ ) can be expressed en- 
tirely via the Gaussian averages of W^s. This 
leads to the following representation of the me- 
ander numbers: 



fc=i 



E 



(4.64) 

Gauss 



where the average over M^'s is calculated with 
the Gaussian weight — the same as in (4.62). 
This formula can be proven by calculating the 
Gaussian integral over (/) in Eq. ( 4.62| ), expanding 
the result in c and comparing with the RHS of 



Eq. (4.63). The factorization at large N is also 
used. 

The principle meander numbers M„ are given 
by Eq. ( 4.64 ) as the linear-in-m-terms, i.e. as lin- 
ear terms of the expansion in m. This looks like 
the replica trick which suppresses higher loops of 
the field W. 

The ordered but cyclic-symmetric sequence of 
indices ai, 02, . . . , a2n-i, Q2n is often called as a 
word constr ucted of m letters. The average on the 
RHS of Eq. ( 4.64 ) is the meaning of a word. Thus, 
the meander problem in equivalent to summing 
over all the words with the Gaussian meaning. 

Since for to = 1 

(2n)! 



\n 



Cm 



(4.65) 



Gauss {n+l)\n\ 
which is known as the Catalan number of the or- 



der n, one gets from Eq. (4.64) 

n 



E- 

fc=i 



(4.66) 



This is nothing but the first sum rule of |2C 



4.2. Complex matrix model for meanders 

The meander numbers can alternatively be rep- 
resented as the Gaussian average over the com- 
plex matrices. The corresponding generating 
function reads 

1 

with 

N 



He) 



In 



E>2 e 



-s 



Gauss 



(4.67) 



S ■ 



-tr 



N ,2 



-cTV^tr (0il¥i02VF,) 



(4.68) 



Here, 0i and 02 are Hermitean while Wa {a = 
1, . . . , to) are general complex matrices. 

It is convenient to introduce one more generat- 
ing function 



M(c) 



t>2 e 



I' 



02 e 



-s 



) (4.69) 

Gauss 



where only one component of Wa, say the first 
one, enters the averaging expression. Differenti- 
ating the generating function (4.67) with respect 
to c and noting that all to components of Wa are 
on equal footing, we get the relation 



dc 



mM{c) 



(4.70) 



between the two generating functions. 

In order to show how the complex matrix model 
recovers the meander numbers, let us replace 



or (f>2 in the numerator of Eq. (4.69) by 



N~^d/d(j)-'i or N^^d/d(j>2 , respectively, and in- 
tegrate by parts. Repeating this procedure itera- 
tively, we get 



Mic 
with 

n 



n=l 



m 



k-1 



(4.71) 



k=l 



E 



fc=i 



1,0271 = 1 



■ Wn VF^ 

a2n-l a2r. 



Gauss 



(4.72) 
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Equation ( 4.72 ) can alternatively be derived by 
calculating the Gaussian integrals over and 02 
in Eq. (4.69) by virtue of 



,e^^ = det-i/2 



a, 6=1 



For TO = 1 the formula 



—tv{ww'^Y 



Cm 



Gauss 



(4.73) 



(4.74) 



which is analogous to Eq. ( 4.65| ), holds for 
the complex matrices. This results again in 
Eq. (g^). 

It is instructive to consider also the case when 
W is a fermionic Grassmann valued matrix d la 
||23|| .0 We shall denote the fermionic matrix as F 
and its conjugate as F. Then we get p3| 



Gauss 

n = 2p (even) 
Cp n^2p+l (odd) 



(4.75) 



Since each loop of the fermionic field is accom- 
panied by a factor of (—1), we arrive at the sum 
rule 



fe=i 



n = 2p (even) 
Cl n^2p+l (odd) 



(4.76) 



This is nothing but the second sum rule of |2C[| . 

Note that the trace of the square of a fermionic 
matrix vanishes because of the anticommutation 
relation imposed on the components. This is why 
we did not consider Hermitean fermionic matri- 
ces and used first a representation of meanders in 
terms of complex matrices to discuss fermionic 
representation of meanders. Fermionic matrix 
models are a natural representation of the notion 
of the signature of arch configurations of pO| . 



5 See 



for a review. 



4.3. Supersymmetric matrix model for 
principle meander 

Having the representation ( 4.71 ), ( [4.72| ) of 



meanders via general complex matrices (either 
bosonic or fermionic), we can utilize the idea of 
supersymmetry to kill the loops of the H^-field in- 
stead of the replica trick. Let us consider the two- 
component Wa whose first component is bosonic 
while the second one is a fermionic matrix: 



Wa^iB,F), Wa = W^^{B^,F) 



(4.77) 



Then all the multi-component meanders in 



Eqs. ( [4.7lD or ( [4.72[ ) vanish and we get the fol- 
lowing representation for the principle meander 

2 

Af„= 

a2 ,■■ ^'3'2ti - 1 ,0-271 —1 

I Gauss 

y~i\MWa,^Wa,^_,---Wa,B\ (4.78) 

\ / Gauss 

where we kept trace of the order of matrices of 
how it appears from Eq. (4.71). The signs are 



essential for fermionic components. 

The generating function ( |4.67 ) equals zero for 
the supersymmetric model since all the loops 
of the B and F fields are mutually cancelled. 
One should use alternatively the generating func- 
tion ( 4.69| ) which can be represented for the su- 
persymmetric matrix model as 



M(c) 



N 



-trBBnn 



t)id(p2 e 



-s 



Gauss 



(4.79) 



where S is explicitly given by 



S = ytr0^ + ytr</>2 - cNtr {(j)iB^2B) 
-ciVtr (0iF02F) 



(4.80) 



as is prescribed by Eq. ( 4.68 ) with Wa substi- 
tuted acco rding to Eq . ( 4.77 ). The equivalence of 
Eqs. ( 4.69| ) and (4.79) in the supersymmetric case 
can be proven replacing B in the integrand on the 
RHS of Eq. (|7|) by N'^d/dB^, integrating by 
parts, and recalling that J-{c) = 0. 
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Equation ( 4.78 ) is a nice representation of the 
principle meander which looks more natural than 
the one based on the replica trick. A hope is that 
it will be simpler to solve the m = 2 supersym- 
metric model than a pure bosonic one at arbitrary 
m. 

The total number of nonvanishing terms on the 
RHS ofEq. (|4J1) for the pure bosonic case, which 
we shall denote as #„, is given by the following 
generating function 



n=0 



2n 



^^l_4(m-l)c2-f + 1 



(4.81) 



This formula can be derived |19| using non- 
commutative free random variables. 

4.4. Relation to the KM model 



Equation ( 4.81 ) is known from the solution 
of the KM model with the Gaussian potential. 
There is the following reason for that. Suppose 
that the matrices Wa are unitary instead of the 
general complex ones. Then one has 



(^trC/„,{/,V..[/,,„_,C/t^„^ 



1 for closed loops 
for open loops 



Haar 



(4.82) 



where the average over the unitary matrices [/'s 
is over the Haar measure. 

Here the loops represent the sequences of in- 
dices {fli, 02, ... , a2n-i, o,2n}- The nonvanishing 
result is only when the loop is closed and encloses 
a surface of the vanishing minimal area, i. e. each 
link of the loop is passed at least twice. This is 
a reflection of the so-called local confinement in 
the KM model. 



The generating function ( 4.81 ) coincides with 
the following correlator in the KM model with the 
Gaussian potential on an infinite D-dimensional 
lattice 



E 

n=0 



2n 



tr02(O) 



(4.83) 



where the average is defined with the same weight 
as in Eq. (3/7), provided that m — 2D. 

The solution of the KM model with the Gaus- 
sian potential can be completely reformulated as 



a combinatorial problem of summing over all such 
closed loops of a given length with all possible 
backtrackings (or foldin gs) in cluded. Its solu- 
tion ||5| is given by Eq. ( [4.8l[ ). 

By virtue of the E guchi-Kawai reduction [2^^, 
the correlator ( 4.83 ) in the KM model on the in- 
finite lattice is equivalent to that in the reduced 
model given by 



n=0 



/d<^id02n:=i'^f/ae-^['^-^l^tr<^2 



with TO 2D and the reduced action being 



(4.84) 



S [0, U] = f tr 



cAT^tr ((/)iC/t02C/a) 



a=l 



(4.85) 



The representation (4.84), (4.85) can now be 
rewritten as 



71=1 



a=l 

X / /'^'^i'^'^2e-^[^-^l^tr0i^[/t02t/a \ 



/d0id(/.2e-^^['>'C^l 



Haar 



since the determinant ( 4.73| ) is equal to a con- 
stant when Wa are unitary. The representa- 
tion ( 4.86D loo ks very similar to the generating 
function (4.69) of the meander numbers. The dif- 
ference is that the average is over the unitary ma- 



trices in Eq. (4.86) and over the Gaussian com- 
plex matrices in Eq. ( 4.69| ). 

We can interpolate between the two cases by 
modifying the weight for averaging over W^s 
along the line of ASM]. Let us introduce 



a I ^2a 



F [W, VF^] . (4.87) 



Then the averaging over the Gaussian complex 
matrices is reproduced as a — > while the average 

^See ]27| for a review. 
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over the unitary matrices is recovered as a — > oo 
since the matrix Wa is forced to be unitary as 
a — > oo. 

We see, thus, that the words are the same both 
for the meander problem and for the KM model. 
The only difference resides in the meaning of non- 
vanishing words — it is equal to one for the uni- 
tary matrices. 

5. Conclusions 

• The Kazakov-Migdal-Penner model is an 
explicit example of crumpled strings (— dis- 
cretized random surfaces) in the D > 1 em- 
bedding space. It possesses the only con- 
tinuum limits associated with lower dimen- 
sional theories of c = or c = 1 . 

• The meander problem results in a more 
complicated matrix model than those 
solved before. It belongs to the same 
generic class of problems of words as the 
large- iV QCD in Z? = 4 but is presumably 
simpler. 

• The relation of the matrix models describ- 
ing meanders with the KM model might be 
a hint on how to solve the former ones. 

• The supersymmetric matrix models of the 
type discussed in connection with the mean- 
der problem could be useful for discretiza- 
tion of super-Riemann surfaces and super- 
strings. 
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